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ABSTRACT 

Five ultrahigh energy cosmic rays in the combined AGASA and HiRes stereo 
data are analyzed to test whether they come from a single source. The quad 
above 37 EeV in the 94-event high energy dataset can be analyzed without con- 
sidering magnetic dispersion. The probability that it is a chance association 
is 10~^. Assuming the source is continuous, the random magnetic deflections 
these UHECRs accumulated en route from their source can be used to estimate 
^J {B'^\)D 1 nG-Mpc. A quintuplet including a HiRes event between 10 and 
30 EeV is well fit with this value. Galactic magnetic deflection appears to be 
smaller than in some models. 

Subject headings: cosmic rays, magnetic fields, cosmology: large-scale structure 
of universe 



1. Introduction 

The combined ultrahigh energy cosmic ray data from the AGASA and stereo HiRes 
detectors was recently studied using the maximum likelihood method by the HiRes collabo- 
ration and this author (Abbasi et al. 2004b). It was noted that a 37.6 EeV HiRes event and 
the previously-reported triplet from AGASA (Hayashida et al. 2000) could have a common 
origin: the distribution of their arrival directions is consistent with a pointlike source given 
the measurement errors of the events. The frequency in random simulations of an equally 
strong clustering signal being observed by chance in a dataset of this size (57 AGASA events 
above 40 EeV and 37 HiRes events above 30 EeV) was reported to be about 0.5% using the 
likelihood ratio alone. In this Letter I point out that a more informative criterion for the 
strength of the cluster gives a chance probability of 10~^, and that the uncertainty or bias in 
the estimate due to definition of the dataset is minor. Thus the evidence that these events 
come from a single source is quite strong, and justifies investigating the implications of the 
cluster assuming it is real. Here, I concentrate on inferences which can be made about the 
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magnetic field in the direction of the source. I also investigate a nearby event in the 10-30 
EeV HiRes data (Abbasi et al. 2004a). It gives an excellent fit to coming from the same 
source, but given the size of the low energy dataset has about a 1 in 6 chance of being a 
random association. 

Empirically, the four higher energy events have arrival directions so close that they can 
be analyzed without including magnetic deflection (Abbasi et al. 2004b). Initially without 
magnetic deflection, I explain the method and estimate the chance probability for this cluster. 
Next, magnetic deflection is added to the analysis and the observed arrival directions arc fit 
to the hypothesis that they are spread not only by measurement error but also by dispersion 
due to random magnetic fields between source and detector. Possible net deflection due to 
the regular local galactic magnetic field is also considered. 



2. Mctximum Likelihood Method 

If the magnetic dispersion along the line-of-sight to a source of UHECRs is small and 
the source is strong enough, the maximum likelihood method (ML) gives an important 
improvement over previous methods of studying clustering. It incorporates the resolution of 
individual events and therefore avoids arbitrariness in the definition of a cluster, it allows 
data from different experiments to be combined, and it is the best way to estimate the 
direction of the source. The application of the maximum likelihood method to this problem 
is described in (Abbasi et al. 2004b). One postulates the existence of a point source located 

— * 

at 9s contributing Ug events to the dataset which has a total of N cosmic ray events, the 
remaining N — Ug events being random, isotropically distributed background events. If the 
ith event is a source event, then its true arrival direction is and the probability density for 
observing it in some direction 9 is denoted by Qi{9, 9s), which depends on the angular errors 

— * — * 

of the event and the angular displacement between 9s and 9. On the other hand, if the ith 

— * 

event is a background event, then the probability density for finding it at a location 9 is given 

— * 

by Ri{9), the relative exposure of the detector to an isotropic background of cosmic rays. 

— * 

Each of these functions is normalized over all directions 9 in the sky. Assuming Ri is small, 

— * 

N is large, and N Us, the partial probability distribution of a given arrival direction 9 for 
the ith event is simply 

m = ^Q^iO,es) + ^^j/^R^{&) , (1) 

and the hkelihood for the dataset of A'^ events is C{ns) = YliLi^ii^i)- The best estimate 

— * 

for the source position 9s is the value which maximizes JC. It is convenient to work with the 
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logarithm of the (dimensionless) likelihood ratio, InTZ: 

lnnn,)^ln^ = y lnh(9l^-l]+l\ , (2) 
^(0) \n[ R^^ei) J i ' 

where £(0) is the hkelihood function for Ug — 0. Below, InTZ and Ug denote the values of these 
quantities at the location and Ug value which maximizes InTZ, unless otherwise specified. 

In the simplifying case of Gaussian resolution with a single variance a, uniform exposure 
density R, average inter-event spacing large compared to a, and N ^ Ug ^ 1, one can recast 
(2) in the form 

InTZ = --x^ - lnPe,n,,N (3) 
where = Ylgource minimized with respect to source location, and 

Pe,n.,N = ^ — (4) 

ng\{N-ng)\ \ngj 

is the chance probability to get Ug hits in N tries, in a region whose fractional size is 
e = 27ra^R. 

Note that £ is a probability density and not a probability, and TZ is not a ratio of 
probabilities. To obtain the ratio of probabilities that the given observation is a signal 
versus being background, one must specify precisely what condition to impose defining the 
signal. For instance, one could define the signal probability to be the fraction of times a 
source of Ug events has a > the one observed, and the background probability as the 
fraction of times a random distribution has Ug events with ^ that observed. Or, one 

— # — * — * 

could integrate C{9) over a region of 9 around the maximum 9g, subject to some condition 
on the desired signal strength. 

Directly computing the signal-to-chance probability ratio by integrating C{ng) and £(0) 
is cumbersome and an alternative approach is commonly used in experimental analyses in 
which many synthetic datasets are created (here, by isotropizing the observed events weighted 
according to the exposure of the detectors), and the number of cases having InTZ > IriTZohs 
is tabulated. However this does not in general provide an accurate estimate of the ratio of 
signal to background probabilities, because, as can be understood from equations (3) and 
(4), since eA^ is small, InTZ can be large either because ng is large, or because is small. 
The mean value of x^ for a genuine source is equal to the number of degrees of freedom, 
and a low value of x^ is not an indicator of a "better" source. Thus a better estimate of 
the chance probability of a given observed cluster is obtained by counting the number of 
cases in random trials in which both InTZ > InTZohs and ng > n^*. For this purpose, the 
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simplest choice is Ug* — ns,ohs, and that is used below. Another option which Auger might 
wish to adopt is to choose (by simulation in advance of the analysis) values of InTZ^, and n^* 
that maximize rejection for purely random datasets while keeping, say, 95% of true clusters 
embedded in random background. This procedure gives a minimum value of n^* which can 
be probed with a given density dataset, for a specified degree of background rejection. 



3. Data and Simulations 

For AGASA, the relative exposure is quite uniform, with R = 0.18 sr~^ = 5.5 x 
10~^deg^^. The HiRes relative exposure is similar but less uniform, as shown in Fig. 1 
of Abbasi et al. (2004b). In the region of the quad, Rhr = 0.2 sr^^. The resolution of stereo 
HiRcs events is, with few exceptions, well-approximated as a symmetric, 2-d Gaussian with 
O'ER = 0.4° (Abbasi et al. 2004b). A definitive study of the clustering of the data must 
await publication of the full AGASA catalog, including individual angular errors, but we can 
make a preliminary analysis, as follows. Individual AGASA events can be described by a 
symmetric 2-d Gaussian (M. Teshima, private communication) whose resolution varies with 
energy and depends on whether the event is well-contained or not. A higher energy event 
can be reconstructed better because more detectors are hit and the statistical fluctuations 
within the detectors are lower. Only detectors near the shower core are used in the angu- 
lar reconstruction, so events whose core is at least 1 km inside the array ("well-contained" 
events) are better measured than those at the perimeter (merely "inside"). 60% of the events 
above 40 EeV are "well-contained" and we denote their resolution by (Jg{E), while 40% are 
merely "inside" with resolution denoted ab{E). Fig. 1 of Takeda et al. (1999), shows the 
opening angles containing 68% and 90% of events as a function of energy up to 100 EeV, 
allowing one to infer (last digit not significant): 

ag{E) = 2.13°- 0.80° logio(^Ecv) (5) 
a,{E) = 4.49° - 1.56° logio(^Eev) ■ (6) 

The three events in the AGASA triplet are all well-contained (M. Teshima, private commu- 
nication), so the analysis below uses ag for those events, i.e., 0.62°, 0.74° and 0.75° for the 
77.6, 55, and 53.5 EeV events respectively. In simulations, errors on other AGASA events 
are assigned at random between (7g{E) and ab{E) in the 60% - 40% proportion and all HiRes 
events have cr = 0.4°. 

There are 251 HiRes events in the fiducial region, —10° < dec < 80°; their arrival 
directions can be taken from the skymap in Abbasi et al. (2004a). In addition to the 37.6 
EeV event appearing in the dataset above 30 EeV, one lower energy event is near the quad. 
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at {171.7° , 57.8°}. Its energy will be denoted E^. Using an E"'^-'^ spectrum, the average 
energy of events in the 10-30 EeV range is 15.4 EeV, so for definiteness £^5 = 15 EeV is 
adopted below. Where relevant, the lack of knowledge of E^ is taken into account. Note 
that the energies of several events reported as above 40 EeV in Hayashida et al. (1996) were 
revised downward; here, as in Abbasi et al. (2004b), the specifics of the AGASA events are 
taken from Hayashida et al. (2000). 

There are three technical differences between the simulations here and in Abbasi et al. 
(2004b). First, to find the most likely source location, Abbasi et al. (2004b) scanned over 
a grid of points with about 0.1° separation and took the position with the largest value of 
InTZ. Here, the best source location is roughly determined by a grid-based scan, then the 
exact point in Ug, RA, and dec which maximizes InTZ is found. Second, Abbasi et al. (2004b) 
took the AGASA resolution to be a sum of two Gaussians, while here individual events 
are given Gaussian errors of different sizes, to more accurately reflect the actual conditions 
in AGASA. The value of InTZ for the actual 94 event dataset with the error assignments 
used here is 12.88, whereas InTZ = 12.98 with the non-gaussian errors used in (Abbasi et al. 
2004b). Third, Abbasi et al. (2004b) include the variation of the HiRes exposure density 
over the sky, whereas here the HiRes exposure is fixed to a constant, 0.2 sr~^, which is 1.07 
times the average unit-normalized exposure. This value is chosen to give the correct rate of 
random 4-fold coincidences, based on approximating the HiRes exposure shown in Fig. 1 of 
(Abbasi et al. 2004b) and using the analytic discussion in the previous section to see that 
the rate of n-fold coincidences ~ / {NagRag{(^) + NhrRhr{())Y<P'()- When HiRes publishes 
its exposure in a more accessible form, it will be possible to improve the precision of this 
estimate. 

Now we turn to the data, and consider different approaches to assessing the likehhood 
that the 4 events above 37 EeV come from a common source. The left and right panels 
of Fig. 1 show results from about 17,000 random datasets composed of the 94 events in 
the combined AGASA and stereo HiRes data above 40 and 30 EeV respectively. With no 
condition imposed on InTZ^ most trials have = 1. Removing these events by requiring 
InTZ > 6 (left panel) shows that almost all of the remaining trials have rig < 3; by requiring 
InTZ > InTZobs the background is dramatically reduced as shown in the right panel. The 
distribution in for a genuine source is very different from the random case: almost all 
the trials have ng greater than ritrue — 0.2. This is illustrated by the center panel of Fig. 1 
which shows the ng distribution for 10^ simulated datasets with an embedded quad: 4 events 
distributed according to the 2-d Gaussians of the quad events, and 90 events distributed 
at random. The cut > obs = 3.89 removes very little of the signal but removes more 
than half the remaining random background. The only condition imposed by Abbasi et al. 
(2004b) was InTZ > InTZohs- Requiring instead that both InTZ > InTZohs and ng > ns,obs, the 
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Fig. 1. — Histograms of rig values as described in the text. 



chance probabihty of the observed quad in about 17,000 datasets is found to actually be 
1.0 ± 0.2 X 10~^, where the error is statistical. 

The chance probability is remarkably insensitive to the angular resolution of the events. 
To illustrate, consider the magnetic dispersion analysis of the next section where a gaus- 
sian random variance 32.5°EeV/i? is added in quadrature to the measurement resolution. 
Repeating the simulations above with the new larger errors gives indistinguishable results 
to within the statistical accuracy used. While this may seem counter-intuitive, the analytic 
discussion of the previous section shows it is to be expected at some level of accuracy: dif- 
ferentiating equation (3) with respect to a, one finds that for the relevant range of A^, n^, a 
and R, the change in cancels most of the change in the background probability term. Of 
course, some changes in the relative errors of different events can make a significant change 
in the likelihood, so this insensitivity to resolution is generic buit not universal. It can also 
be understood by noting that the chance probability of finding 4 random hits in a clustered 
configuration depends only on how the spacing of the clustered hits compares to the average 
density of events, and has no relation whatever to the accuracy of the measuring apparatus, 
which is solely relevant for determining the quality of the fit to the source hypothesis. As 
a result of this natural insensitivity to the resolution, the conclusion that the quad events 
have a common source is unlikely to change when more precise resolutions become available 
for individual AGASA events. 

The chance probability of 0.1% is not strictly an a priori probability because the choice 
of the energy threshold of the original AGASA dataset used by Hayashida et al. (1996) was 
not specified before that analysis, and the correct relative normalization of energies between 
AGASA and HiRes is uncertain. However these are minor issues in determining the chance 
probability of the quad. Since the energy of the events is not inherently important to the 
cluster analysis, the threshold is only relevant to the density of background events. The 
original AGASA doublet events which are members of the quad have energies 77.6 and 55 
EeV - well above the 40 EeV threshold for the AGASA dataset, and the HiRes member 
with energy 37.6 EeV is well above the threshold of 30 EeV so there is no issue that the 
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background is being unnaturally minimized by these thresholds. (Analyses relating to the 
total number of doublets or higher multiplets are sensitive to the definition of the dataset, 
since one of the three original AGASA doublets had a 43 EeV event so that shifting the 
threshold say to 50 EeV significantly changes the conclusions (Finley & Westerhoff 2004).) 

There is another measure of the probability that the quad is real, which is completely 

independent of the dataset choice. One can ask how often starting with the initial AGASA 
dataset of 30 events (Hayashida et al. 1996) and adding 27 more AGASA and 37 HiRes 
events at random, does one find InTZ > InTZohs, "^s > "^s, obs- The result is 32 cases in 5 x 10'^ 
tries, for a chance "promotion probability" of 0.7 ± 0.1 x 10^'^, of which more than half of 
this comes from the "CI" AGASA doublet, which was not in fact promoted. 

There are many "interesting" anomalies which could occur in a random dataset, and the 
question comes up whether to sum over the probability of getting all possible combinations of 
712 doublets, ns triplets, etc, when computing the random background, which would increase 
the background estimate. The correct procedure depends on the theoretical perspective. 
If all sources are of comparable apparent luminosity and magnetic smearing is negligible, 
multiple clusters are expected and their presence or absence can be used to validate or 
exclude such a model. By contrast, for the picture inspiring this study - that different 
sources probably have substantially different apparent luminosities and cosmic magnetic 
fields are probably very inhomogeneous, so that few if any sources can be reconstructed - 
one has no a priori expectation to find multiple significant sources and the confidence in the 
hypothesis is not increased by having multiple sources. Rather, the presence or absence of 
multiple reconstructed sources gives information on the nature of the sources and fields. In 
this scenario, the important question is how well the observed cluster fits the hypothesis; 
the probability that the cluster in question could happen by chance determines how much 
credence should be attached to its implications while awaiting further data. 



4. Magnetic Deflection and Dispersion 

The strength and direction of the local galactic magnetic field is known (Heiles 1996) 
but its decay with distance from the local spiral arm and the partition between random and 
coherent components are quite uncertain. The measured local field may well be dominated 
by the random component, yet its direction might be correctly modeled, since global informa- 
tion constrains the direction of the coherent field for a given assumed structure. Therefore, I 
consider three cases for the galactic magnetic field below: a) Taking the coherent magnetic 
dcficction to be negligible, appropriate if the local field is dominated by the random com- 
ponent, b) Taking the direction of galactic magnetic deflection to be flxed as in standard 
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models but the magnitude of the galactic deflection, G*/E^ to be free. This is appropriate 
if the global constraints correctly determine the general direction of the coherent field along 
the trajectory but, for instance, not how rapidly it decreases away from the disk, c) Taking 
the direction of deflection as well as to be free. 

In the limit of sufficiently many small deflections, the smearing of cosmic rays from a 
pointlike continuous source due to random fields is described by a 2-d Gaussian probability 



where A is the coherence length of the field. D is the distance to the source, and (B^X) denotes 
a weighted mean of B^X over the trajectory (Waxman & Miralda-Escude 1996). I assume 
throughout that all events in the cluster are protons; there are too few events to permit 
generalizing the analysis to different charges at this time. Since the experimental angular 
resolution of each event is also approximately a symmetric 2-d Gaussian, of width denoted 
(7i,o, magnetic smearing increases the effective size of the individual events' resolutions in 



with the same charge and energy follow the same trajectory and have no relative dispersion. 
Since a caustic only operates over a very narrow range of energy it is a priori improbable for 
any given energy, so the relative magnetic dispersion between the 55 and 53.5 EeV events 
should be minimal. Therefore when determining the magnetic smearing we first combine 
these two events, replacing them by a single event with arrival direction {170.22°, 56.57°}, 
energy E = 54.25 EeV, and a = 0.52°. 

It turns out that cases a) and b) are practically the same, because when the direction of 
the coherent deflection is fixed toward —38.6° in local equatorial coordinates, both quad and 
quint are fit best with = 0. This is also true for case c) when all five events are used; it 
is not meaningful to study this case with just the quad, since there are 6 measured positions 
and 5 parameters to be constrained, so the uncertainty is too large. Given the inherent 
uncertainties, G* ~ is consistent with known constraints on the strength and extent of 
the local coherent field, as will be developed in greater detail elsewhere. The magnitude of 
deflection predicted by standard models gives a poor flt to the observed arrival directions of 
the quad, and a very poor flt for the quint (G. Medina- Tanco, private communication). 

Now the magnetic smearing can be estimated. The simplest procedure is to adjust 
to the value which causes for the cluster to be equal to the number of degrees of freedom 
in the flt, as expected on average. One can also choose the value of E^ which maximizes 
InTZ; in general this gives a lower value of E^ since smaller a means lower background. We 
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Barring a caustic, two cosmic rays from the same source 
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adopt the former because it is computationally easier but keep in mind that the true value 
of may be lower. With = 0, there are three free parameters: the coordinates of the 
source and E^,. The "quad" consisting of the 77.6, "54.25", and 37.6 EeV events is therefore 
described by a distribution for 6 - 3 = 3 degrees of freedom, and a first estimate of E^, 
is that which gives = 3 for the "quad", namely E^, = 32.5°EeV. Fortunately, E^, is quite 
insensitive to a possible systematic offset in the energies between AGASA and HiRes events. 
Each experiment has an energy uncertainty of ~ 15% and a commonly evoked explanation 
for the discrepant spectra is that AGASA and HiRes energies are systematically 15% too 
high and low respectively. (Indeed, awareness of this possibihty is why the HiRes threshold 
of 30 EeV was investigated in (Abbasi et al. 2004b).) But rescahng the AGASA energies by 
0.85 and the HiRes energy by 1.15, leads to E^ = 32.3°EeV. 

If the fifth event does indeed come from the same source, including all 5 events in the 
analysis will give a more sensitive determination of both because the number of events 
increases so the distribution becomes narrower and because the low energy event is more 
sensitive to E^. However to help judge whether the 5th event is likely to be from the same 
source, we first compute fo^^ the "quint" using the value E^ = 32.5°EeV obtained from 
fitting the quad. The result, = 3.5, is a better-than-average fit since there are now 5 
degrees of freedom. The fact that ioi the quint is so good, using the E^, fit to the quad, 
is evidence that the 5 events are really from the same source and that our treatment of 
magnetic effects is robust. If the fifth event had instead been a random background event it 
most probably would have fallen as far as possible from the quad source location, since that 
is where the phase space is largest, and the E^ required to get a reasonable fit to the quint 
would have been larger than found with just the quad. 

Estimating E^ from the "quint" as the value for which x^ — 5, gives E^ — 23.9°EeV and 
a source location of {169.40°, 56.78°}, with a 90% error radius of ^ 1°. This is the same 
as that obtained with the quad, {169.24°, 56.81°}, within the uncertainty. Even though 
the 5th event gives a better-than-average fit to the point-somcc hypothesis, decreasing the 
HiRes energy threshold to 10 EeV increases the dataset from 94 to 308 events, so there 
is a significant (1 in 6) chance the fifth event is an accidental association. The promotion 
probability from doublet to quint is 0.16% - small but twice that from doublet to quad - 
because in this case 251 instead of 37 HiRes events are added to the original 30 and added 
27 AGASA events; the promotion probabihty is insensitive to the value of E^ taken here to 
be 23.9°EeV. 

Assuming the 5 events have a common source, we obtain a "90% CL" range for E^ 
as follows. For 5 degrees of freedom (continuing to merge the 53.5 and 55 EeV events), 
10% of the time > 9.2 and 10% it is < 1.6. The values of E^ giving those values of x^ 
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are — 12°EeV and 55°EeV respectively. The analogous range using just the "quad" is 
(5 — 96)°EeV, confirming the stronger constraint which is possible when all 5 events are used. 
If the energy of the fifth event is lower (higher) than the 15 EeV used above, its effective a 
would be greater (less) for a given E^,, changing the value of i?* derived. Taking the energy to 
be 10 EcV, the best value and "90% CL" ranges become 20(9 - 52)°EeV, while for E^ = 30 
EeV, E^ = 30(18 — 65)°EeV. Until the energy of the fifth event is made available, we can 
best characterize E^ as 23.9°EeV with range (9 — 65)°EeV, always remembering that the 
maximum likelihood results are somewhat lower. 

Equation (7), which translates into E^, = \J2{B'^\)D /?>, is valid only when the net 
deflection angle times distance to the source is large compared to the typical coherence 
length of the field: E^/E x D > X. There is a strong deficit of L* or brighter galaxies in the 
direction of the UHECR cluster, out to 140 Mpc (G. R. Farrar, A. A. Berhnd and D. W. 
Hogg, 2005, in preparation). The range of coherence lengths for extragalactic fields usually 
seen in the literature is ~ 0.1 — 1 Mpc, so with A = Ampc Mpc and D = 140 D^q Mpc, 
the condition of validity is < 70 /^mo/Ampc EeV. Our deductions about the magnetic 
field structure should therefore be trustworthy if the source is beyond the void, particularly 
since lower energy events are most crucial in determining E^,. However if the source is 
from one of the unremarkable, low- luminosity galaxies around D ^ 26 Mpc, our results 
are only qualitatively valid unless A < 0.2 Mpc. Note that an additional criterion must 
also be met if the Gaussian-random-variable approximation is to be accurate. Namely, the 
various members of the cluster apart from the merged 53.5 and 55 EeV pair, must sample 
independent and uncorrelated regions of magnetic field en route to Earth. This requires 
£"^1-^ — -^1 > \/D, which is more stringent than the condition on the highest energy event 
alone but may be satisfied at least marginally. As more events become available, it will be 
necessary to improve upon the analysis of Waxman & Miralda-Escude (1996) and include 
correlations between the defiections of events of neighboring energies, unless X/D is very 
small. 

Keeping in mind the caveats above, we convert our best value E^, = 23.9°EcV to an 
estimate of the integrated extragalactic fields. (It is qualitatively clear, given the very small 
scales involved in the random galactic fields compared to the Larmour radius of a proton, 
Rl — 108^^^ kpc, that deflection by random galactic flelds is comparatively insigniflcant 
as was recently shown in Tinyakov & Tkachev (2004).) From eqn (7): 

^{B^X)D < 0.96 (0.36 - 2.6) [1.3] nG Mpc , (8) 

where the uncertainty due to the energy of the 5th event not being available is given in 
parentheses and the value is square brackets is that derived from the quad alone. If the 
source is in fact in a galaxy beyond the void which extends to 140 Mpc in the direction of 
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the cluster, we can obtain an estimated upper limit on the rms field strength in the foreground 
void: (-Bvoid) ^ 0.07/-y/AMpc nG. This is compatible with Faraday Rotation estimates from 
distant quasars, as shown by the sheets- voids model of Farrar & Piran (2000). It is also 
below the equipartition value, very crudely estimated by taking the void material to be fully 
ionized, with an electronic density 10% of the mean baryonic density, and a temperature of 
T3 lO^K as suggested by Furlanetto et al. (2004), giving Egqui = 0.3 T^^^ nG. 

5. Summary and Conclusions 

Four cosmic ray events in the combined stereo HiRes and AGASA data above 30 and 40 
EcV arc consistent with having a common origin, with a chance probability of 10~^. When 
the HiRcs threshold is lowered to 10 EeV, the datasct increases from 94 to 308 events and 
a fifth event consistent with coming from the same source is found. However due to the 
large number of added events, there is a 1 in 6 chance that it is background. The best 
estimate of the source direction is {169.3°, 56.8°} in J2000 equatorial coordinates, with 
an 1° 90% error radius. No coherent net defiection, expected from galactic magnetic 
fields, is required. The integrated magnetic smearing is used to infer the order of magnitude 
constraint a/ {B'^X)D ?a 1 nG Mpc, valid for sufficiently large D or small A. If the source is 
beyond the void which extends to a distance of 140 Mpc (Farrar et al, 2005, in preparation), 
we obtain an upper limit on the magnetic field in the void: (-Bvoid) ^ 0.07/ y^Xupc nG. With 
the low statistics available at this time, these constraints on the intervening magnetic fields 
must be considered qualitative rather than quantitative, but it is interesting that magnetic 
fields have such small effects in at least one arrival direction. 

The existence of a single source of 4 or 5 events excludes models in which UHECRs are 
products of the decay or annihilation of superheavy relic particles, unless a mechanism can 
be found to cause an extreme and concentrated overdensity of rehcs along this line-of-sight 
which seems highly unlikely. This cluster also adds to the difficulties of the Z-burst model, 
because within that model such a cluster could only reasonably occur if produced by the 
same (tremendously powerful) source of neutrinos, presumably at a cosmological distance. 
But in that case, there should be negligible magnetic smearing within the cluster since the 
observed UHECRs would be mainly photons from 7r°'s produced in a local neutrino DM 
halo. 

A probability of ~ 0.1% that the quad is due to a statistical fiuctuation corresponds 
to its being a "3.3 cr" effect. This is well above threshold for being taken seriously but is 
not definitive. More data will be very helpful - higher energy events to establish the signal 
more firmly (or debunk it) and lower energy events to permit a more detailed study of the 
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source and magnetic field. HiRes is continuing to run and AGASA plans to release its lower 
energy data soon, however these data will make a fairly meager change in the statistics. The 
tremendous physics payoff of having much better statistics on this cluster provides additional 
motivation for a new northern hemisphere UHECR detector with large aperture and good 
resolution since it is not in the field of view of the Pierre Auger Southern Observatory. 
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